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ABSTRACT 

A mechanism of generation of a toroidal component of large scale magnetic field, 
leading to the reconstruction of the pulsar magnetospheres is presented. In order to 
understand twisting of magnetic field lines, we investigate kinematics of a plasma 
stream rotating in the pulsar magnetosphere. Studying an exact set of equations de- 
scribing the behavior of relativistic plasma flows, the increment of the curvature drift 
instability is derived, and estimated for Is pulsars. It is shown that a new parametric 
mechanism is very efficient and can explain rotation energy pumping in the pulsar 
magnetospheres. 
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1 INTRODUCTION 

The aim of the present work is to investigate generation of 
a toroidal component of the magnetic field nearby the light 
cyhnder surface (LCS) (a hypothetical surface, where the 
linear velocity of rotation equals the speed of light). 

The work we consider in this paper is closely related to 
the pulsar wind problem. Studying the magnetic field of the 
Crab nebula, Piddington (1953) was first who has suggested 
the presence of a central object in the nebula, with frozen 
magnetic field inside. It has been supposed that rotation of 
the central body provokes generation of the toroidal com- 
ponent of magnetic field. Further investigations have shown 
that this kind of magnetic field characterizes magnetized 
star winds (IWeber fc Davis 1967|) . These results have been 
generalized for relativistic flows in a region close to the LCS: 
UMichel 19691 IKennel et al. 19831 IKennel fc Coroniti 19841 
[Bejel man fc Li 1992p . Despite success of developed models, 
they encounter a number of difficulties, when one attempts 
to extrapolate the wind back to the source: the pulsar mag- 
netosphere. For large distances the wind is specified in the 
approximation: a = B'^ /(iirmc^n'y) <^ 1, where B is the 
magnetic field induction, m and n- the electron mass and 
density respectively and 7 the Lorentz factor of relativis- 
tic electrons. In this case, change of magnetic field's con- 
figuration is defined only by plasma motion. This circum- 
stance simplifies a possibility of analytical consideration of 
a plasma. But in the pulsar magnetospheres, a situation is 
opposite, the energy density of magnetic field exceeds by 
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many orders of magnitude the energy density of the plasma 
(T S> 1, therefore a need of consideration of this specific case 
is essential. Close to the light cylinder area the magnetic field 
drags behind itself the rotating electron-positron plasma and 
the question which arises is: how the magnetosphere is re- 
constructed nearby the light cylinder surface? It is obvious 
that close to this region, rigid rotation is impossible and 
consequently magnetic field lines must deviate, lagging be- 
hind the rotation of the pulsar. Implementing special MHD 
codes in a series of works UMichel fc Krause-Polstorff 19841 
IKrause-PolstorflF fc Michel 19851 ISmith et al. 2001 P pulsar 
wind physics has been numerically studied and improved 
by Spitkovski fc Arons (2002) and Spitkovski (2003) where 
plasma dynamics in 3D was presented and it has been shown 
that the flow goes through the LCS into the wind zone. In 
these papers a principal assumption is the current generated 
by the electric drift: Vi = x 'S/B'^ (Blandfo rd 2002) . 
Obviously for a plasma composed of equal numbers of posi- 
tive and negative charges, the current is not generated (the 
electric drift does not " feel" charges) , although for the pulsar 
plasma a primary electron beam is composed of only elec- 
trons and therefore the electric drift generates the current, 
leading to creation of electromagnetic fields. 

In (iRo gava et al. 2003[ ) a particle moving along a curved 
rotating channel has been considered and it was shown that 
for a certain shape of curved trajectories one may avoid 
the light cylinder problem. Therefore one has to understand 
what is a mechanism responsible for the process of twisting 
of field lines when the condition ct S> 1 is satisfied. 

According to observations it is clear that the energy of 
emission is very high. An observed pulsar luminosity lies in 
the range: [10^^ - 10^^]erg/s (|Tores fc Nuza 2002 |l . on the 
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other hand the only source of pulsar radiation can be ro- 
tational energy Ilu'^/2, where I is moment of inertia of the 
pulsar, and uj - the angular velocity of rotation. As observa- 
tions show the spin down luminosity is of the same order of 
magnitude as the radiation luminosity, therefore it is reason- 
able to suppose that all pulsars emit due to rotation energy 
decrease l|Sturrock 1970|) . The problem concerns the ques- 
tion: how the rotation energy is transformed into pulsar ra- 
diation. According to standard models, due to electric field, 
the charged particles uproot from a surface of the neutron 
star and accelerate by the electric force which results in the 
radiation process. The origin of this emission is supposed to 
be in the magnetosphere of pulsars. These models introduce 
a vacuum gap, inside of which the particles experience strong 
electric field and accelerate. But the problem arises concern- 
ing the gap size which turns out to be not enough for energy 
gain of charged particles (|Ruderman fc Sutherland~1975 , ) . 

In order to resolve this problem and enlarge the 
gap size (which will provide increase of an acceler- 
ation length scale) many attempts have been done, 
applying different approaches: fArons fc Sharleman 1979| 
[Mushmov fc Tsygan 1992| IRuderman fc Sutherland 1 9751), 
but no approach was able to get the efficient acceleration 
enough for producing observed radiation. 

A new mechanism of acceleration has been introduced 
in ([Machabeli & Roga va 1994p where a bead moving inside a 
straight rigidly rotating pipe has been studied. It was shown 
that the centrifugal force can be very efficient and if one ap- 
plies this method for the pulsar magnetospheres it will pro- 
vide high Lorentz factors of particles. Therefore the amount 
of energy contained within the e+e~ plasma is very high. If 
one finds mechanisms for the conversion of at least a small 
fraction of this energy into the variety of waves or instabil- 
ities - one might witness a number of well- pronounced and 
bona fide observational signatures in the pulsar radiation 
theory. In ( [Machabeli fc Rogava 1994[ ) it has been found 
that the radial component of velocity for relativistic particles 
behaves in time as ccos{^lt) (c is the speed of light), which 
gives a possibility of parametric energy pumping from the 
mean fiow into instabilities (see (Mac habeli et al. 2005|l ). In 
(|Machabeli et al. 2005P the e"'"e~ plasma has been studied 
and the increment of an instability of the Lengmuire waves 
was estimated. It has been demonstrated that the centrifu- 
gal acceleration might have been efficient enough for the 
observed spin down luminosity. We have shown that the lin- 
ear stage was so efficient that it was very short in time, and 
nonlinearities were turned in soon. 

In the present paper we generalize the previous work 
and study the parametric mechanism of the curvature drift 
instability driven by the centrifugal acceleration. We con- 
sider a two component plasma: a) the basic plasma flow 
(bulk flow) with the concentration Upi and the Lorentz fac- 
tor 7p; and b) the beam component with the concentration 
n(, and the Lorentz factor 71,. It is known that in the pul- 
sar magnetosphere the drift velocity is to be important for 
plasma dynamics. The drift velocity may influence processes 
in the plasma and especially may affect an evolution of insta- 
bUities. Unlike ( .Spitkovsky fc Arons 2002||Spitkovsky 2003P 
where the processes are considered nearby the pulsar sur- 
face, in the present paper we investigate instabilities close 
to the light cylinder area, where effects of centrifugal accel- 
eration should be extremely efficient. In ( [Spitkovsky 2003| 




Figure 1. Here we show geometry in which we consider our sys- 
tem of equations. By e^, er and ex unit vectors are denoted, note 
that ex -L e^.^. C is the curvature center. 

it has been noted that the structure of pulsar magneto- 
spheres could not be solved analytically, whereas in the 
present paper, we show that an initial stage of the recon- 
struction process of magnetospheres can be considered ana- 
lytically, starting by appropriate initial conditions. Another 
difference is that in our model we study a plasma, which 
is bound by rigidly rotating straight magnetic field and the 
force free condition applied in ( [Spitkovsky fc Arons 2002] 
[Spitkovsky 2003 [) is not valid, because as it is shown in 
( Shapakidze et al. 2000[ ) the force free condition can be pro- 
vided only if the magnetic field has a configuration simi- 
lar to the one of a differentially rotating Couette ffow. The 
principally different assumption in the present paper is that 
instead of considering the electric drift, we study the cur- 
vature drift investigating the possibility of generation of the 
toroidal component Br, which is a key step in understanding 
the reconstruction of the pulsar magnetosphere nearby the 
LCS. 

The work is organized as follows. In !j2]we derive the dis- 
persion relation, in Sj31the corresponding results are present 
and in ^we summarize the results. 



2 THEORY 

Throughout the work it is supposed, that magnetic field lines 
are almost straight and due to the frozen-in condition the 
plasma particles follow the magnetic lines and accelerate. 
Geometry in which we consider the problem is shown in 
Fig. d] ( [Lyutikov et al. 1999| ) . 

Our system is governed by the Euler equation 
l|Machabeh et al. 2005| ): 

+ (viV)pi = -7«Va + - (E + Vi X B) , (1) 
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(14) 



i = pl,b, 



A^{t) = — sin(m), 



the continuity equation: 



dt 



V(n,Vi) = 0, 



(2) 



and the induction equation, which closes the system: 

VxB = if + ^5:j, (3) 
c at c ^-^ 

where J, (i = pi, b) is the current of plasma and beam com- 
ponents. 

We start our analysis by introducing small deviations 
around the equilibrium state: 

+ (4) 

where * = (n, v, p, E, B). 

Since we are interested in the generation of the toroidal 
component of magnetic field, it is interesting to study the 
curvature drift wave (when lj ~ ku^), because it is char- 
acterized by the following conditions: Br 3> -B^, ^ Er 
( [Kazbegi et al. 1991[ ) that show an importance of Br on the 
one hand and closes the system by the second condition on 
the other hand. 

We consider the equilibrium state with a drift velocity 
along the a;-axis 

U)Bt(.B 

where 7 is the Lorentz factor, ujb = e_Bo/(m-ec) and 
Rb is the curvature radius of the magnetic field lines 
(e and m are the charge and mass of electron and Bq- 
the magnetic induction). Along (j), due to the centrifugal 
acceleration one has a relativistic flow with the velocity 
( |Machabeh fc Rogava 1994[ ): 

uo0 — ccos(r2t). (6) 

If one expresses the perturbation of physical quantities 
by following: 



*^(f,r) (X 'I'^(t)exp[j (kr) 



(7) 



then considering only x components of the Euler and induc- 
tion equation, it is easy to show that for curvature drift 
waves, propagating perpendicular to magnetic filed lines 
(k^ <^ kx), Eqs. pi2l3|l can be reduced into the form: 



dPix -r, , \ 1 6 7->l 
t(kxUOx + K^Uo^)Pix — -V04,Br, 

— i{kxUox + kcj,uo4,)nl — ikxTiiovl, 



ik^cBl = 47re > [rnovlx + nlviOx). 



E 

i—pl.b 



(8) 

(9) 

(10) 



In Eq. ^ we have used an approximate expression of 
velocity along the r-axis: Vr ~ cEx/Bq^. If we choose p\x 
and n] to have the form: 



where 

A4t) = ^ + ^sin(2m), 



(11) 
(12) 

(13) 



LObRb ' 

then one obtains from Eqs. (|8I9|I 



(15) 



vL = / e-*^''*')^;o4^')B.(^')d^^ (16) 



^1 _ «enoifc^^ikAi(t) / 



'VQ^{t )Br(t )dt . 



(17) 

Substituting Eqs. (|16ll7p into Eq. (fTOjl it reduces to the 



form: 



k^cBl{t)= y ^e*Ai(*) f e-^'^^^^''^vo4t')Br{t')dt'+ 
,K 7iO J 

2 rt rt" 

— kxUoixe ' / dt / 



i—pl,b 



vo<p{t )Br(t )dt , 



(18) 



where uji — ^ Anmoe^ /m is the plasma frequency. In or- 
der to solve this equation one has to take the Fourier time 
transform. For this reason if one uses the following identity: 



itia; sin Qt \ ^ t 

e -2^Js 

s 

one can reduce Eq. (fTSi 



BriLo) ^ y y y 

i—pl,h cT — il s,n,l,p 



xBr {lj + Q (2[s -l] + n-p + a)) 



(19) 



Ja{gi)Ju{h)Ji{gi)Jp{h) 
u+i^^ +n{2s + n)' 



Uix 



E 



t^i kxUix 



E E 



Js(gi)Jn{li)Ji(gi)Jp{h) 

xBr{u) + Q,{2[s-l + fi\ + n-p + a)), (20) 
where 



gi = 



kx U 

k^c 



3 DISCUSSION 

One can see from the dispersion relation that the system 
is characterized by two different kinds of resonance, which 
come from the first and second terms of the right hand side 
of Eq. ([20)1: 



(21) 
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Figure 2. Dependence of the increment on A nearby the LCS. 
The set of parameters is P = Is, 75 ~ 10*, = lOOORjc- 



and 

w+:^^+n(2[s + M]+n) ~0, 



s,n = {0,±l,±2...}, 



(22) 



±1. 



As we will see later, a resonance frequency from the first 
condition Eq. (|21|) does not influence the corresponding res- 
onance from the second expression, because if the first term 
is valid, the second condition is not satisfied and vice versa, 
when the second resonance works, the first one is not valid. 

Since we are studying the mechanism of energy pumping 
from pulsar's rotation and a natural consequence, that mag- 
netic field lines must be twisted nearby the light cylinder and 
the twist should have a direction opposite to the rotation, 
one can suppose that a frequency responsible for this process 
must be small compared the angular velocity of the pulsar. 
On the other hand, assuming A ~ 6 x 10^" — 3 x 10^^ cm 
(where A « Aa, is the wave length) for second pulsars it is 
straightforward to check that \kxUix/2\ ~ 0.02 — 0.1, which 
for any values of s and n is much less than i}{2s + n) for 
s,n 7^ (see Eq |2Hl . Thus the only possibility which pro- 
vides low frequency waves from the first resonance condition 
is: 



2s -f n = 0. 



(23) 



Here we have assumed that Rb ^ Rp ^ 10^ cm {Rp is 
the pulsar radius). Unlike this case, second resonance condi- 
tion does not provide low frequencies (see Ea. (l22|l '). because 
even for vanishing s and n, jj, is not vanishing and hence it 
does not contribute in the process of magnetic field line's 
twisting. 

Let us consider the dispersion relation near the beam 
resonant condition expressed by Eq. (I21|l . Then only reso- 
nant terms will be preserved and Eq. p0|l will reduce: 



LOi,kxUbx 



EE 



Br[uJo) 

CT — it 1 S,l,p 

X Br {iiJo + ^ {-21 - p + O-)) , 

where 

kx Uhx 

ijJo ~ t: — 



Js(gb)J-2s{h)Ji{gi)Jp{h) 



A2 



(24) 
(25) 



and the frequency has been expressed by the form: 

Lo = ujo + iA. (26) 

Here A's imaginary part A = /m(A) is related to the 
increment of the instability. Since a dominant term in Eq. 
(|24[) comes from low frequencies {ujo <^ Q), then the only 
terms contributed in a time average will have p equal to 
Q{2{s — l)+n + a), because all other terms with Br{uJo + ^q) 
{q 7^ 0) give zero due to an oscillative character with very 
big values of frequencies. Taking into account this condition, 
one gets: 



73/ \ _ ^bkxU, 

Br (UJO) 



■.Ubx \ ^ \ ^ 

kj,c ^ ^ 



Js{gb)J~-23{h)Jl{gb)J-2l + a{h) 



bOl^cp 



o- = ±l s,l 



A2 



X Br{iUo) 



(27) 



From here one can easily express the increment by fol- 
lowing: 



'^'''-^'"'E^(gb,h)E2{gb,h) 



2'ybok^c 
where 

Ei{gb,h) = ^ Js{gb)J-2a{h) 

S 

^2{gb,h) = J2 ^Jiigb)J- 



(28) 



(29) 



(30) 



Strictly speaking Ei and E2 are functions of gb and h 
and one can show, that these summations are convergent 
(see Appendix). 

It is interesting to investigate the increment versus fol- 
lowing physical quantities: Xx (for the fixed, and very big 
in comparison with Xx values of A^). On the other hand 
one has to compare results with an observational evidence. 
As we have already mentioned the only source that may 
provide energy for radiation is the slowdown of the pulsar: 
W « IQCI, here / is moment of inertia of the pulsar. The 
rate of rotation energy loss can be estimated by following 
ratio: W /W ~ 2Cl/n = 2P/P, where P is rotation period 
of the pulsar. The given ratio is different for different pul- 
sars and ranges from 10~" sec"^ (PSR 0531) to 10"^* sec"^ 
(PSR 1952-f29). Therefore the increment of the instability 
must not be less than 2P/P. 

We investigate the instability rate nearby the light 
cylinder, because the centrifugal acceleration should be most 
efficient in this region. In Fig. [5] we show dependence of the 
increment on: A^; nearby the LCS. The set of parameters is 
P — Is, jb ~ 10*, Xx ^ X, X^ — lOOOiijc and it is supposed 
that kx < and Ubx > (otherwise the resonance frequency 
is unphysical - negative). Here Ric is the light cylinder ra- 
dius. Such a choise of X^ provides almost perpendicular (to 
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the equatorial plane) propagation of waves. One can see that 
the increment reaches the value 10~^, which is more by 
many orders of magnitude than typical values of 2P/ P. The 
linear stage will be very efficient and short in time strongly 
indicating the non linear regime of the phenomenon. A need 
of non linear saturation is seen also from the fact that Br, 
which is responsible for twisting of magnetic field lines, oscil- 
lates with frequency ljo, due to this oscillation, Br not only 
will lag behind the rotation, which is physically reasonable, 
but also will advance it, therefore the need of non linear sat- 
uration of the instability increment is essential. Therefore 
initially created small perturbations will rapidly increase in 
time and thanks to the instability process, it will extract 
energy from the background flow into the energy of electro- 
static waves. 



4 SUMMARY 

(i) Considering the relativistic plasma flow composed of 
the primary and secondary (beam) components, we have 
studied the role of the centrifugal acceleration in the curva- 
ture drift instability. 

(ii) Making the linear analysis of equations, we have de- 
rived the dispersion relation and a new mechanism of the 
parametric instability responsible for rotation energy pump- 
ing has been found. 

(iii) Considering low frequencies, which are responsible 
for twisting of magnetic field lines, an expression for the 
instability increment has been obtained. 

(iv) Studying dependence of increment on Xx, it has been 
found that the instability was very efficient and increments 
were more than pulsar spin down rates by many order of 
magnitude indicating the need of the non linear considera- 
tion of the problem. 

As we have seen, the analysis indicated the importance of 
the non linear stage in dynamics of the instability, therefore 
it is essential to study the same problem numerically by 
implementing a special relativistic MHD code, which will 
comprise one more step closer to the real scenario. 
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APPENDIX A: 

In this section we would like to show that the sum: 
Ei(ff6,/i) = ^J,(3i,)J„2,(/i), (Al) 

s 

is finite. 

In order to prove the convergence of (A.l) we use the 
following inequality: 



i^^(^)|^t>Vt)(I)^^^"'' 

Here: v > — \- 

Then for the term in T,-i[gb, h), one can write: 



(A2) 



2s+l 



(A3) 



(s + l)!(2 

where we have used the well known equivalence: 

J_,(x) = {~lYJ,{x). (A4) 



This condition shows that Cs = 
where: 



Us = 



(s-f l)!(2s + l) 



(f)(l 



\Js(gb)J2s{h)\ < Us, 



(A5) 



Let us prove that the Us is convergent using the Dalam- 
ber criterion, by introducing the following ratio: 



Us 



u 



^ - (f ) (i 



(s + 2)(2s + 2)(2s + 3)' 



It is obvious that one can find so for which 



(A6) 



= q< 



1. We see that the expression in Eq. (|A6|l decreases with an 
increasing value of s, i.e: 

2 1 



U 



Us 



< q - 



(f)(l 



(so + 2)(2so-h2)(2so-|-3) 



< 1, (A7) 



S > So. 



This means that the sequence Us is convergent, and 
hence for s > the summation Ei((/6,/i) is flnite. 

When considering the case s < and formally intro- 
ducing a new index m = — s, one obtains an expression: 



\Jm{gb)J2m,{h)\, 



(A8) 



similar to a corresponding term in Eq. (|A3|l for s > and 
hence, the summation for negative values of s is also con- 
vergent. 

The proof for convergence of the second summation 
T,2{gb,h) (see Eq. (|3U|l 'l does not principally differ from the 
one we have already considered and therefore we do not show 
it here. 
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